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Abstract 

This paper is concerned with the asymptotic behavior of solutions of 
stochastic differential equations dyt = duit — W{yt)dt, yo = 0. When 
d = 1 and V is not periodic but obtained as a superposition of an infi- 
nite number of periodic potentials with geometrically increasing periods 
{V{x) — 'Y^'^=oUk{x/ Rk), where Uk are smooth functions of period 1, 
J7fe(0) = 0, and Rk grows exponentially fast with k) we can show that 
yt has an anomalous slow behavior and we obtain quantitative estimates 
on the anomaly using and developing the tools of homogenization. Point- 
wise estimates are based on a new analytical inequality for sub-harmonic 
functions. When d > 1 and V is periodic, quantitative estimates are ob- 
tained on the heat kernel of yt, showing the rate at which homogenization 
takes place. The latter result proves Davies's conjecture and is based on 
a quantitative estimate for the Laplace transform of martingales that can 
be used to obtain similar results for periodic elliptic generators 



1 Introduction 

It is now well known that natural Brownian Motions on various disordered 
or complex structures are anomalously slow. 

These mechanisms of the slow diffusion for instance are well understood 
for very regular strictly self-similar fractals. The archetypical specific ex- 
ample of a deep problem being the one solved in l |Barlow and Bass, 1999| l 
on the Sierpinski Carpet (which is infinitely ramified, a codeword for hard 
to understand rigorously: for a survey on diffusions on fractals we refer to 
HBarlow, 1998 1 , for an alternat ive approach to jOsada, 1995^ and for the 
random Sierpinski Carpet to ( jHambly et al., 1998[ ). It appears that the 
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main feature is the existence of an infinite number of scales of obstacle 
(with proper size) for the diffusion. 

It is our object to show that one can implement the common idea that 
this last feature (infinitely many scales) is the key to the possibility of 
anomalous diffusion, in a general context using the tools of homogeniza- 
tion. 

The strategy of the proof might appear paradoxical: it is not a priori 
very sensible to try to prove that the diffusion is anomalous by the use 
of homogenization theory which is a vast mathematical machine destined 
to prove an opposite result, i.e a central limit theorem and thus normal 
diffusion. But it will be shown that when the homogenization process is 
not finished, an anomalous behavior whose characteristics are controlled 
by homogenization theory might appear. 

This paper will focus on the sub-diffusive behavior in dimension one (sub- 
section which will allow the introduction of a concept of differentia- 
tion between spatial scales that can be applied to a more general frame- 
work. 

The proof of the anomaly of the exit times is based on a new quantitative 
analytical inequality for sub-harmonic functions (subsection 12.311 that is 
linked with stability properties of elliptic divergence form operators. 
The extension of those results to higher dimensions has been done in 

( |Ben Arous and Owhadi, 2001| l and to the super-diffusive case in | |Ben Arous and Owhadi, 2002^ 
and HOwhadi, 2001 b> . 

The control of the anomalous heat kernel tail is based on sharp quantita- 
tive estimates for the Laplace transform of a martingale. These estimates 
allow us to put into evidence the rate at which homogenization takes place 
on the behavior of the heat kernel of an elliptic generator in any dimension 
f subsection 12.211 . The quantitative control of the heat kernel in homog- 
enization theory outside any asymptotic regime has been recognized as 
difficult and important | |Norris, 1997^ . For instance, this problem is at 
the center of Davies' conjecture emphasized as "well beyond existing re- 
sults" (|n avies, 19931 . With theorem l2.8l we give a proof of that conjecture 
in any dimension for elliptic operators with only bounded coefficients. 



1.1 History 

The idea of associating homogenization (or renormalization) on large num- 
ber of scales with the anomaly of a physical system has already been ap- 
plied from an heuristic point of view to several physical models. 
Maybe one of the oldest of such applications is to Differential Effective 
Medium theories which was first proposed by Bruggeman to calculate the 
conductivity of a two-component composite structure formed by succes- 
sive substitutions (( [Bruggerman, 1935| l and I |A1P, 1977} ) and generalized 
m orris, 1985^ to materials with more than two phases. For instance this 
theory has been applied to compute the anomalous electrical and acoustic 
properties of fluid-saturated sedimentary rocks {Sen et al., 1981) . More 
recently this problem has been analyzed from a rigorous point of view in 
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The heuristic application of this idea to prove the anomalous behavior of 



2 



diffusion seems to have been done only for the super-diffusive case that 
is to say for a diffusion evolving among a large number of divergence-free 
drifts. Maybe this is explained by the strong motivation to explore con- 
vective transports in turbulent flows which are known to be characterized 
by a large number of scales of eddies. The first observation was empirical: 
in ( [Richardson, 1926^ Richardson empirically conjectured that the diffu- 
sion coefficient Dx in turbulent air depends on the scale length A of the 
measurement. More recently physicists and mathematicians have started 
to investigate on the super-diffusive phenomenon (from both heuristic 
and rigorous points of view) using the tools of homogenization or renor- 
malization (the first cousin of multi-scale homogenization): we refer to 
(Avellaneda and Majda, 199 0'^; fGlimm and Zhang , 1*9921 , jAvellaned a, 1996|l; 



( Bhattacharya, T999f7lFWnjiang and Papanicolaou7^994ll; fFannjiang^nd Komorowski, 2001||. 



1.2 The model 

Let us consider in dimension one a Brownian motion with a drift given by 
the gradient of a potential V , i.e. the solution of the stochastic differential 
equation: 

dyt = dujt-\7V{yt)dt, yo = 0. (1) 

The multi-scale potential V is given by a sum of infinitely many periodic 
functions with (geometrically) increasing periods: 

oo 

^ = E (2) 

n— 

In this formula we have two important ingredients: the potentials Uk 
and the scale parameters Rk- We will now describe the hypothesis we 
make on these two items of our model. 

1. Hypotheses on the potentials Uk 
We will assume that 

Uk G C°°(T) (3) 

f7fc(0) = (4) 

Here C°°(T'*) denotes the space of smooth functions on the torus 
T :— R/Z. We will also assume that the first derivate of the Uk are 
uniformly bounded, i.e. 

Ki ■— supsup |(7fc(a;) - Uk{y)\/\x - y\. < oo (5) 

k€N x^y 

We will also need the notation 

Ko~supOsc{Uk) (6) 

fceM 

where the oscillation of Uk is given by Osc((7) := sup U — inf U. 
We write D{Uk) for the effective diffusivities associated to the po- 
tentials Uk' if zt is the solution of dzt — dut — 'VUk{zt)dt it is well 



3 



known fOlla, 1994|l that as e J. 0, tz^i^i converges in law towards a 
Brownian Motion with covariance matrix D{Uk) given by 

D{U„) = ( / e^^"(^'da; / e-^^"(^)da;)-\ (7) 
Jt Jt 

We also assume that the effective diffusivity matrices of the Uk's are 
uniformly bounded away from and 1. 

Amin = inf D{U„) > and A^ax = sup D{U„) < 1. (8) 

2. Hypotheses on the scale parameters Rt 

Rk is a spatial scale parameter growing exponentially fast with k, 
more precisely we will assume that Rq — vq — 1 and that the ratios 
between scales defined by (we write N* the set of integers different 
from 0) 

rfc = Rk/Rk-i G N* (9) 

for fc > 1, are integers uniformly bounded away from 1 and oo: we 
will denote by 

Pmin := inf rfc and pmax ■■= sup (10) 

fcSN* fcgfj. 

and assume that 

Pmin > 2 and pmax < OO. (11) 

Since || VV^Hoo < oo it is well known that the solution of Q exists; is 
unique up to sets of measure with respect to the Wiener measure and 
is a strong Markov continuous Feller process. 

Remark 1.1. Note that if Wn,U„ G {Wi, . . . ,Wp} , the (Wi) being non 
constant, then the conditions and are trivially satisfied. 

2 Main results 

2.1 Sub-diffusive behavior 

Our first objective is to show that the solution of Q is abnormally slow 
and the asymptotic sub-diffusivity will be characterized in three ways: 

• as an anomalous behavior of the expectation of r(0, r) (the exit time 
from a ball of radius r, for r oo, i.e. Eo[r(0, r)] ^ r^^"). 

• as an anomalous behavior of the variance at time t, i.e. Eo[i/?] ~ t^^'^ 
as f — > oo. 

• as an anomalous (non-Gaussian) behavior of the tail of the transition 
probability of the process. 

More precisely there exists a constant po{Ko, Ki, Amax) such that 
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Theorem 2.1. // pmin > po o,nd t(0, r) is the exit time associated to the 
solution of Q) then there exists a constant Ci depending on Ko,Ki such 
that 

Eo[r(0,r)] =r^+''i<'''+^<''' (12) 
where e(r) as r ^ oo and 

0<-- <VY{r)<~- h . (13) 

in Pmax Pmin in Pmax in Pmin Pmin in Pmin 

Theorem 2.2. // pmin > po o,nd yt is a solution of Q) then for there 
exists a constant C2 depending on Kq,K\ and a time to depending on 

-K"!, Pmin, Pmax, Amax SUch that for t > to 

nvt] = i'-^ (14) 

where 

InAmax C2 , , InAniin C2 , , 

0<-l 1 i <Mi)<--, r^. (15) 

mpmax In Pmin In Pmax m Pmin (mpmin)"' 

Theorem 2.3. // pmin > po and yt is a solution of mil then there exist 
constants C5 depending on Ko,Ki,R2, C3 on Ko, Ki, pmin, d^Ca^Cj on 
Ko,Ki such that ift,h>0 and 

1>C, and y > C3(-^)^'"^'— ('""-n'^ (16) 

then 



with 



\nr[\yt\>h]<-CeYijir (17) 
In A max 

V3 = > 0. (18) 

In Pmax In Pmin m Pmax 

Remark 2.4. The second condition in is really needed since the lead- 
ing exponent associated to {t/h) is [ 2^^^'^^ ), i-e. half the one associated 
to V2,. This condition corresponds to a frontier with a heat kernel diagonal 
regime. 



2.1.1 Description of the proofs 

Before discussing the results further we want to describe the proof. A 
perpetual homogenization process takes place over the infinite number of 
scales 0, . . . , n, . . .. The idea is to distinguish, when one tries to estimate 
II2I 1. 1141 1 or 1171 1. the smaller scales which have already been homogenized 
(0, . . . , Tie/ called effective scales), the bigger scales which have not had a 
visible influence on the diffusion {ridri, ■ ■ ■ ,00 called drift scales because 
they will be replaced by a constant drift in the proof) and some interme- 
diate scales that manifest their particular shapes in the behavior of the 
diffusion (rie/ + 1, . . . ,ndri — 1 = n-ef + Uper called perturbation scales 
because they will enter in the proof as a perturbation of the homogeniza- 
tion process over the smaller scales). To estimate 1121 for instance, if one 
considers the periodic approximation of the potential 

n 

yo"(a;) = 5](7fc(a;/7?fc) (19) 

fc=0 
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the corresponding process y^"' will have an asymptotic (homogenized) 
variance Holla, 1994| l 

^(^o") = ( [ e2^«"(«"^>da; f e-'"'o"(«-.-)rf^) "i (20) 
Jt Jt 

D{Uo) is smaller than 1 and because of the geometric growth of the peri- 
ods Rn and a minimal separation between them (i.e. pmin > po), -D(Vo") 
decreases exponentially fast in n. 

By homogenization theory, y" is characterized by a mixing length ^m(Vo") ~ 
Rn such that if one writes r" its associated exit times then for r > ^m(Vo") 

Eo[r"(0,r)]~-^. (21) 

Writing ne/ (r) = sup{7i : R„ < r} one proves that Eo [t(0, r)] ~ Eo[r"<=^ (0, r)] 
by showing the stability of Eo[t(0, r)] under the influence of V^^(^j^.j^ — 
X^fc^n /(r-)+i Uk{x/Rk)- This control is based on a new analytical inequal- 
ity which shall be described in the sequel and allows to obtain that 

Eo[r"=^«(0,r)]e-'°=^'-'''"^/M+^' 

(22) 

< Eo[T(0,r)] < Eo[r"<=-f('-)(0,r)]e'°"'''''"=/<-'+i\ 

In these inequalities OsCr(V„^j,(r)+i) stands for sup3(Q K'^j,(r)+i-inffl(o,r) K'^^f^j+i 
and is controlled by 

Osc,(l/^^(,)+i) < Osc((7„^^„+i) + ||VV;^^(,)+2||ocr 

i.e. Uefir) -\- 1 acts as a perturbation scale and n^f^r) + 2, ... ,00 as drift 
scales. From this ^ 

Eo[r(0,r)]~ (23) 

Thus, if 

-lim inf -LinD(lz"=^'''') > 

r^oo Inr 

one has sub-diffusivity, in the sense as defined above. 
The proof of I14II follows similar lines by the introduction mixing times 
Tm(V^) and visibility times r„(T/p°°) (such that for TmiVg) < t < Ty{Vj^), 
V^°° has not a real influence on the behavior of the diffusion yt and V^" 
has been homogenized). Then choosing nej{t) = sup{n : TmiVo) < t} 
one obtains that 

Proposition 2.5. Let i'2{t) the function associated to \i4\l , one has for 

^e/(i)(l - r^) < Mt) < '^efim + (24) 

In Pmin m Pmin 



In- 



i^e/(i) = 1^^^^ With p"'/ =R„^j.and X^^^ ^ D{Vo^'). (25) 
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This proposition shows than this separation between scales is more 
than a conceptual tool, it does reflect the underlying phenomenon. In- 
deed the anomalous function i'2{t) is given in the first order in l/(lnpmin) 
by the number of effective scales by E[j/f] ~ tD{VQ °^), and in this ap- 
proximation V2{t) ~ i^ef{t) where Vefit) corresponds to a medium in which 
the ratios r„ and the effectives diffusivities D{U„) have been replaced by 
their geometric mean over the rie/ -I- 1 effective scales. The origin of the 
constant /(In pmin) in 12411 is the perturbation scales. More precisely, 
one has to fix the drift scales by ndri{t) = inf{n : Tv{V^) > t}, and in 
general there is a gap between n^fft^ and ndri{i), the scales Un situated 
in this gap manifest their particular shape in the behavior of V2{t) and 
since no hypothesis have been made on those shapes one has to take into 
account their influence as a perturbation. 

One may notice that in many papers on diffusions on fractals (see e.g. 
| |Barlow, 1998| l section 3) obtaining estimates on hitting times is essen- 
tially the key to the whole problem and the same is true here: this strat- 
egy has been adapted in l |Ben Arous and Owhadi, 2001| l. In this paper we 
have chosen to not use this strategy in order to put an emphasis on the 
role played by the never-ending homogenization process taking place on 
these diffusions on fractals. Indeed one might wonder why the estimates 
of the behavior of Brownian Motions on fractals are of the form 

E[y?]~t^, (26) 



E[r(0,r)]~r''™, (27) 

lnp(t,a:,y)~-( '^~^''"' )^. (28) 

One explanation is given here by the number of effective scales hidden 
in the estimates I26II . 1271 and 1281 . Let us assume the model to be 
self similar (for all k, rt = p and Uk = U, D{Uk) ~ A). In the table 
below we have summarized formulae giving (in the first approximation 
in 1/ In p) the number of effective scales and the formulae linking them 
with those anomalous estimates (appearing in the proof, the infiuence of 
the perturbation scales will be neglected). This gives three values of 
corresponding to 1261 . 1271 1. 1281 1 and the interesting point is to compare 
them. 

Let us observe that the multi-scale homogenization techniques gives 
back the right forms for the mean squared displacement, the exit times 
and the transition probability densities; they are explained by the number 
of scales which homogenization can be considered as complete associated 
to each observation. Moreover dw^i,djn,2 and dw,3 are equal up the first 
order approximation in 1/ In p nevertheless they are not equal and this is 
not surprising. Indeed when p is small the second order term in l/(lnp)^ 
can not be neglected since the perturbation scales becomes more and more 
dominant (and the infiuence of the perturbation scales is of the order of 
l/(lnp)^). 
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MVt] Eo[t(0,2;)] lnPo[2/t>/i] 
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2 
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A3 
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2.1.2 Strong overlap between the spatial scales 

The anomaly is based on a minimal separation between spatial scales i.e. 
Pmin > po and one might wonder what happens below this boundary. The 
answer will be given on a self similar case, i.e. V is said to be self similar 
if for all n, Un = U and p^in = Pmax = P- 

Theorem 2.6. // the potential V in ^j) is self similar. Then for all p > 2 



with 



Eo[r(0,r)] = r2+''W 
In p 



(29) 
(30) 



with e(r) ^0 as r ~* oo. 

Here Vp is the topological pressure associated to the shift operator 
Sp-. x£T-^px£T (see 112911 for its definition). 

Using the convexity properties of the topological pressure one has 'Pp{2U) + 
Vpi-2U) > and 

Proposition 2.7. Vp(2U) + 'Pp{-2U) = if and only if 



lim II iy (U{p''x) 

n^oo II n ^ 



Yd 



U{x)dx) 



= 0. 



(31) 



From this one deduces that for the simple example U{x) = sin(x) — 
sin(81a;), E[r(0,r)] is anomalous (sub-diffusive ~ j-^+i^ with u > 0) for 
p G {2} U {4, . . . , 26} U {28, ... 80} U {82, . . . , +00} and normal (~ r^) for 
p = 3,27,81. 

Thus if U is not a constant function, there exists po{Ko, Ki, D{U)) such 
that for p > Po, yt has a clear anomalous behavior (Eo[r(0, r)] ~ 
with 1/ > 0) but in the interval (l,po] both cases are possible: yt may 
show a normal or an anomalous behavior according to the value of the 
ratio between scales p and the regions of normal behavior (characterized 
by proposition |2|2j might be separated by regions of anomalous behavior. 
What creates this phenomenon is a strong overlap or interaction between 
scales: that is why the region (l,po) will be called "overlapping ratios". 
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i.e. in this region the fluctuation of at a size f > is not represented by 
a single Un{x/R„) but by several ones and to characterize the behavior 
of yt in that region one must introduce additional parameters describing 
the shapes of the fluctuations f/„, elsewhere a normal or a sub-diffusive 
behavior are both possible. 

2.2 Davies's conjecture and quantitative estimates 
on rate of convergence towards the limit process in 
homogenization 

The proof of theorem 12.31 has not been described yet. The strategy is 
still to distinguish effective, perturbation and drift scales nevertheless it 
is not obvious to determine how many scales have been homogenized in 
the estimation of foivt > h). The answer is directly linked with the rate 
at which the transition probability densities associated with a periodic 
elliptic operator do pass from a large deviation behavior to a homogenized 
behavior. 

Consider for instance in any dimension d > 1, [/ G L°=(T'*) and the 
Dirichlet form 

£:(/,/) = 1/2/ iV/(a.)i^ " /G2?[f] = (32) 

Write p{t, X, y) its associated heat kernel with respect to 

the invariant measure associated to 13211 . Note that when U is smooth 
the associated operator can be written L = 1/2A — V(7V and it is well 
known that 

• Large deviation regime: for |x — j/l >> t the paths of the diffusion 
concentrate on the geodesies and 

lnpit,x,y)^-\^^. (34) 

• Heat kernel diagonal regime: for \x — y\'^ « t, the behavior is 
fixed by the diagonal of the heat kernel and 

pit,x,y)^^. (35) 

t2 

Davies conjectured that (we refer to ( |Davies, 1993| l, he considers periodic 
operators of divergence form nevertheless the idea remains unchanged) 
that p{t,x,y) should have a homogenized behavior {lnp{t, x,y) ~ —{x — 
y)D{U)~^{x — y)/{2t)) for t large enough. 

J. R. Norris | |Norris, 1997} has shown that the homogenized behavior of 
the heat kernel p{t, x, y) corresponding to a periodic operator on the torus 
T'' (dimension d side 1) starts at least for tint >> \x — j/p (with \x — 
y\ « t ); in this paper it will be shown that it starts for t » \x — y\ in 
any dimension. 

This allows to complete the picture describing the behavior of pit, x, y) 
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• Homogenization regime: for 1 << \x — y\ << t and \x — y\'^ >> t, 
homogenization takes place and 

lnp{t,x,y) ~ -ja:: - y\%-^^^/{2t) (36) 

with 

la^ - vId-Hu) ■= \x - v)D{Uy^{x - y). (37) 
More precisely we will prove that 

Theorem 2.8. Consider pit, x, y) the heat kernel associated to the Dirich- 
let form 13 with respect to the measure mu ■ Then there exist constants 
C, C2 depending only on d and Osc(f/) such that for 

C\x-y\<t, C^/t<\x-yl C<\x~y\ (38) 

one has 

p{t,x,y) < 1 e^p {-{1- E)\y-x\l-,^^^/i2t)) (39) 

(27rt)#(det(D(U)))2 

p{t,x,y)> 1 rexp(-(l + £)|y-x|^-i(^)/(2t)). (40) 

(27rt)#(det(D(U)))5 

With 

E{t,x,y):^C2{^-^^ + -^)<^. (41) 
t \x — y] 10 

Theorem EH 

proves Davies's conjecture, moreover E{t,x,y) acts as a 
quantitative error term putting into evidence the rate at which homoge- 
nization takes place for the heat kernel, and it also acts as the inverse of 
a distance from the domains associated to the large deviation regime and 
the heat kernel diagonal regime. Observe that all the constants do depend 
only on d and Osc(i7). It is straightforward to extend those estimates to 
any periodic elliptic operator. They can be liked to results obtained by 
A. Dembo (Dcm bo, 1996^ for discrete martingales with bounded jumps 
based on moderate deviations techniques. 

2.2.1 A note on the proof of theorem 12.31 

Those estimates basically say that the homogenized behavior of the heat 
kernel associated to a periodic medium of period R starts for t > R\x — y\. 
Thus in the proof of theorem 12.31 the number of the smaller scales that 
can be considered as homogenized is fixed by n^fit/h) — sup^{7?„ < t/h}, 
which (assume D{Un) = A and R„ — p" for simplification) leads to an 
anomaly of the form 

lnP(y. >h)< ^C^f-jj- ~ -C^(i)-fe^ ~ ^ 

(42) 

with ~ 2 — The equation l|42|l suggests that the origin of the 

anomalous shape of the heat kernel for the refiected Brownian Motion on 
the Sierpinski carpet can be explained by the formula linking the number 
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of effective scales and the ratio t/h. 

The first condition in lltill can be translated into "homogenization has 
started on at least the first scale" and the second one into "the heat kernel 
associated to Q is far from its diagonal regime" (one can have h? /t « 1 
before reaching that regime, this is explained by the slow down of the 
diffusion) . 



2.2.2 A quantitative inequality for exponential martingales 

The core of the proof of theorems 12.31 and 12.81 is an inequality giving a 
quantitative estimates for the Laplace transform of a martingale: 
Consider Mt a continuous square integrable adapted martingale such 
that Mo = and for A G R, E[e^*^*] < oo. Assume that there exists a 
function / : such that for all t2 > > one has a.s. 

t2 t2-tl 

E[J d<M,M >s iJ^tJ < J f{s)ds (43) 

tl 

with f{s) = /i for s < to and /(s) = f-z for s > to with to > and 

0</2 </l. 



Theorem 2.9. Let Mt be the martingale described above. 

1. for all < |A| < (2e(/i - /2)to)"^ one has 

E[exp(AMt)] < e3(i-i/9(^» exp(^AV2t) (44) 

with g{\) = i-x'Hjl-f2)to^ ^^"-^ verifies l<g<2 

2. for alio < ly < (2e(/i — f2)to) one has 

E[exp(;. <M,M >0] < ^M^htf-^^^^P^^^r-^- (45) 

((/i - tv^^oY 

This theorem uses the knowledge on the conditional behavior of the 
quadratic variation of a martingale to upper bound its Laplace transform, 
and it is well known that a quantitative control on the Laplace transform 
leads to a quantitative control on the heat kernel tail. The condition A 
small enough marks the boundary between the large deviation regime and 
the homogenization regime. A direct application of the key theorem is the 
following result. 

Corollary 2.10. Let Mt be the martingale given m theorem \2.iA 

Write Ci = (2e(/i - /2)to) V/a- For r = ^ < I one has 

2 

r{Mt>x)<e^'-\^p{-{l-r^)^). (46) 

2/21 

This corollary gives a quantitative control on the tail of the law of Mt 
from the asymptotic behavior of its conditional brackets. 
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2.3 An analytical inequality for sub-harmonic func- 
tions 

The stability property 1221 is based on the following analytical inequality: 

Theorem 2.11. Let Q. be an open bounded subset of M. (d = 1), for 

X G C°°(n) such that X > on Q, and (f>,ip € C^{fl) null on dfl and both 
sub harmonic with respect to the operator — V(AV), one has 

I X{x)\\/(l>{x).\/tp{x)\dx <3 I X{x)\/(l>{x).\/tp{x)dx. (47) 
Jn Jn 

The constant 3 in this theorem is the optimal one. We believe that 
this inequality might also be true in higher dimensions, i.e.: 

Conjecture 2.12. For Q C R'^an open subset with smooth boundary, 
there exist a constant Cd,n depending only on the dimension of the space 
and the open set such that for X G C°° (SI) such that X> Q onQ, and 4>,ip & 
C^(r2) null on SJl and both sub harmonic with respect to the operator 
~V(AV), one has 

I X{x)\\/(l>{x).\/ip{x)\dx<Cd,n I X{x)\/(l}(x).\/iP{x)dx. (48) 
Jn Jn 

This conjecture is equivalent to the stability of the Green functions 
of divergence form elliptic operators under a deformation. More precisely 
write Gx the Green function associated to — V(AV) with Dirichlet condi- 
tions on dQ. 

Proposition 2.13. The coniecture \2.1'A is true with the constant C'd,n if 
and only if for all A,/i bounded and strictly positive on 

(supmax(^,A))-^^.- < ^^pyl < (supmax(^,^))'''''". (49) 
n A ' Gx{x,y) ^ X fi ' 

Remark 2.14. Thus it would be interesting to prove it since it would 
allow to obtain sharp quantitative estimates on the comparison of elliptic 
operators with non Laplacian principal part. By proposition 12.131 it is 
easy to check that conjecture 12 . 1 21 implies Harnack inequality. One might 
think that one would be able to obtain 14911 using Aronson's estimates and 
keeping track of the dependence of the constants in the Harnack inequality, 
but this is not the case since Harnack inequality is an isotropic inequality 
and 1491 compares in an optimal way Green functions of operators which 
can be strongly anisotropic. 

Let us remind that the Harnack inequality associated to the operator 
L = — VAV says that for all L-harmonic functions u in 5(0, r) one has 

sup u{x) < Cl inf , 

xeB(0,r/2) xeB{0,r/2) 

where the optimal constant Cl grows towards infinity as sup A/ inf A ^ oo 
whereas the constant associated to conjecture 12.121 is independent of A. 
That is why the Harnack inequality strategy, which has already been used 
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to obtain quantitative results for the comparison with the Laplace opera- 
tor ( we refer to ||Sta mpacchia, 1965| l, l |Ancona7 "19971 , jGriiter and Widman, 1982l l 
and | |Pinchover^^89| l) allows to obtain 

Pf^, < Ch (50) 

but with a constant Ch exploding like Cd exp (Cd(sup A/ inf A)*^"*) 
Remark 2.15. Since the conjecture is true in dimension one with Cd,n ~ 3 
(this constant is an homotopy invariant), it is through proposition 12.131 
that one obtains stability property (I22II . 

Remark 2.16. It is easy to deduce from Theorem 12.111 that if f2 is a 
bounded open subset of R"^ and 4>, tp are both convex or both concave 
functions on and null on dfl, then 

/ \V:,(l>{x).V^i^{x)\dx<3 [ V^(l>{x).V^i^{x)dx. (51) 
Jn Jn 

Remark 2.17. The conjecture 12. 121 (theorem 12. Ill when d = 1) has an in- 
teresting signification (and consequences) in the framework of electrostatic 
theory, we refer to the chapter 13 of ( |Owhadi, 2001a[ . 



2.4 Remark: fast separation between scales 

The feature that distinguishes a strong slow behavior from a weak one 
is the rate at which spatial scales do separate. Indeed one can follow 
the proofs given above, changing the condition pmax < oo into Rn = 
Rn^i[p" /Rn^i] {p,OL > 1) and Ainax = A,nin = A < 1 to obtain 

• A weak slow behavior of the exit times 

Cir^e"'") < Eo[T(0,r)] < Car^'e^*"' (52) 

with g{r) = (In r)i (In 1/A)(lnp)" = . 

• A weak slow behavior of the mean squared displacement 

Cite-f^'^ < Eo[yt] < Cafe"^'*' (53) 

with f{t) = (lnt)»(lnl/A)(21np)-»(l-he(t)) 

• A weak slow behavior of the heat kernel tail: for h > 0, Ci < t/h < 
Cafe 

¥[yt >h]< C3e~'^^^'''-T^^ (54) 

with fc(a;) = A-(t^'°<i+'("» 

And as Q I 1 the behavior of the solution of Q pass from weakly anoma- 
lous to strongly anomalous. 
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3 Proofs 



3.1 Davies's conjecture and quantitative estimates 
on rate of convergence towards the limit process in 
homogenization 

3.1.1 Quantitative control of the Laplace transform of a 
martingale: theorem 12.91 

The core of the proof of the anomalous heat kernel tail (theorem 12.311 
and the quantitative estimates on the heat kernel associated to an elliptic 
operator (theorem 12.811 is theorem 12.91 that will be proven in this sub 
subsection. 

Let Mt be the martingale described in theorem 12.91 Let g > 1, using 

2 

Holder inequality and Ito formula it is easy to obtain that with hq — 2(q-i) 

E[exp(AMt)] < ElexpihqX^ <M,M >t)]^ (55) 

Thus the quantitative control of the Laplace transform of the martingale 
shall follow from this control on its bracket. 

Write ^ = ^ ([n] shall stand for the integer part of /i). Using Holder 
inequality and the control 1431 one obtains for 1 < z < oo 



E[exp(/igA^ <M,M >()]? <E[exp(z/i,A^ < M,M >[^]to)]''' 

exp [(hq/q)\^{t - Mto)/i)- 
Then by taking the limit z [ 1, one easily obtains that 

E[exp(/i5A^ <M,M >t)]' <E[exp(/i,A^ < M,M >[^,tj]9 

exp {{hq/q)X^{t- [^]to)fi). 



(56) 



(57) 



Write a= we will need the following lemma 

Lemma 3.1. Let Mt he the martingale described in theorem \2.fA and rj > 

0, for a — /2//1 and fi = t/to one has 

E[exp(r7<M,M>0]<l + 2^^^4r^ E it^ - m)" Cl^ (a - ir . 

n—l 0<7n<nAfi 

(58) 

Proof. By the Taylor expansion of the exponential one obtains 

+00 

exp(?? <M,M>t) = l + Y, V"W„ (59) 

n=l 

with Wn = J 1{0 < ti <■■■< t„ < t)d < M, M ---d < M,M >t„. 
Using the control I43II on the conditional brackets of the martingale it is 
easy to obtain by induction on the integrand and the Markov property 
that 

E[W„] < J l(0<uiH \-u„ <t)f{ui)--- f{u„)dui---du„ 
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Combining this with I59II and using the fact that /(s) < fig(s/to) with 
g{z) = l(z < 1) + al{z > 1) one obtains that 

+00 

E[exp(ry < M, M >t)] < 1 + ^ (r?/ito)"G„ (60) 



with G„ = / 1(0 < 21 + ■ ■ ■ + 2„ < At) n (l(2fc < 1) + al{zk > 

l))dei ■ ■ ■ dzn- Developing the product in Gn one obtains by integration, 
induction and straightforward combinatorial computation that 



□ 



G„ = -!- V C™(Ai-m)"(a- l)"". 
n\ ^ — ' 

Which leads to (I58II by the inequality 16011 . 

Using lemma one obtains 

E[e:^p{hgX^ < M,M >Mto)] 

^glMlfitoT J2 (M-m)"Gr(a-l)'" 

n=0 ' 0<m<nA[ii] 

Changing the order of summation, one obtains 

E[exp(/i,A^ < M,M >[^]t„)] < exp(^AVitoM) 

y- (M-m)'"(fe,(a-l)AVito)'" (61) 
^ m! 

0<m<[fi] 

Now we will need the following lemma 
Lemma 3.2. for — i < y < 

(M-m)"'y"' ^ exp(yM) 

0<m<[M] 

Proof. Put — i < a:: < and write for n G N, 

Ex™ 
— r(n — m)™ 
m! 

0<m<ri 

It will be shown here that Vp G N* , Vn G N 

/ri < (itp(a;)) " (1 — itp exp(a;Up)) ^ (63) 

where ttp the increasing sequence defined by = and Up+i = exp{—xup) 
and converging to j/o the smallest positive solution of yexp{xy) — 1. 
The inequality I62II is then obtained for Up{y) = 112(2/) = exp(— y) and 
using exp(— y) — 1 > —y and — i < y < 0. 
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Write y\ = inf {y > : yexp(|a;|i/) — l} (note that < j/i < 1) and 
consider for —y\ < y < yi the function 

/ : y ^ (1 -yexp{xy)y^ 

By Taylor expansion, for y G /(y) = J2n=o T.t=o ^"m^ 

and since 

y^n^ y"-'^^^^^T^ — = -; "'",11. < oo with a normal convergence 

of the series, the order of the limits can be changed, which leads to 

+00 , +00 +00 ^ +00 

f(y) = E ^ E = E ^ E -^rv- 

m—O n— m— n — m 



\-oc n 



En \ ^ / \m \ ^ n r 



n— m—O 



It follows that Vn G N, /„ = -^'"'/"^ Now, for < a; < 0; the 
constant yo = inf {y > : yexp{xy) — 1} does exist and Vy G 
] - yi,yo[, Vn, /("'(y) > (thus 7„ > 0). 

Thus from the classical theorem of Taylor expansion, the series 

^n=o ^ converges towards f for y g] — yi, yo[ and in that interval 

00 

E y"^"^ = {^-V exp(a;y)) ~ 

n — O 

From which one deduces that Vy G]0, yo[ Vn G N In < y "(l— yexp(a;y)) ^ 
On the other hand if one considers the sequence uo = 0, = exp(— a;itp) 

then it is an exercise to show that Up is increasing and will converge to- 
wards yo, which leads to 1631 . □ 

Applying I62II to I6II I with y = hq{a — f\to one obtains that for 

< |A| < {ehq{fi - f2)to)~^ one has 

E[exp(ftqA'' <M,M>[^jtJ]i < exp(^AV2toM) (^(1 - a)AViio)"i 

^ (64) 
Writing v — X^hq and combining 1641 with 15711 one obtains the inequality 
(HSl of theorem IT^ 

Combining 1641 with 1571 and 1551 one obtains the inequality 1441 of 
theorem l2.9l bv choosing q — (A^(/i — /2)toe) ^ (? > 2 under the condition 
imposed on A). 



3.1.2 Upper bound estimate (|^ of theorem EH 

Theorem can be used to give quantitative estimates on any operator 
as soon as a cell problem is well defined. Consider yt is a diffusion on R"* 
that may be decomposed for t > as 

yt=x + x{t) + Mt (65) 



16 



where x(i) is a uniformly (in t) bounded random vector process (Hxlloo < 
Cx) and Mt is a continuous square integrable J-'t adapted martingale such 
that Mo = 0. 

Assume that for all / £ R** with |Z| = 1 there exists a function / : 
such that for all t2 > ti > one has a.s. 

E[J d< M.l, M.l >, \ J^tA < J f{s)ds (66) 

tl 

With /(s) = fi for s <to and /(s) = HDI < fi for s > to with to > 0. 
where D is a positive definite symmetric matrix. 

Assume that the diffusion yt has symmetric Markovian probability densi- 
ties p{t, X, y) with respect to the measure m{dy) such that for all x,y €'R'' 
and t > 

Pit,x,y)<^ (67) 

and for 5 > 

'^o:{\yt-x\>5)<C3e-'^*^ (68) 
where C2,C3,C4 are constants. 



Theorem 3.3. Let yt he the diffusion described above. 

Thenmthki = 30(e(/i-A,nin(-D))to) VVin(-D) andfe = 30+10dAmax(D)(l+ 
ki\x-y\<t, fca < \x-y\>AC^ (69) 

one has 

p{t,x,y)<^e.v[-{l-E) ^y-f^-' ) (70) 

with E, = C2(5(A,,in(D)C4)-^ + 2'*C3) anrf £ = 3((^;i^)2 + < 
10 

Proof. The estimate on the heat kernel p{t, x, y) will follow from the chain 
rule and decomposing it the probability of moving away from x to "a 
well chosen set containing y in the time tg" and its complement. More 
precisely, writing ey-^ :~ {y — x)/\y — x\ and As — {z £ : [z — 
x).ey-x > (1 — 5)\x ~ y\}, using 1671 one obtains that for t > 0, a;, 1/ G R"* 
and < (? < 1, 



P(i, x,y)= / p{tq, X, z)p{t{l -q),z, y)m{dz) 
J As 

+ P{tq,x,z)p{t{l- q),z,y)m{dz) 



C2 r 1 

'---T -j'^AVtq-e^y-^ > I2; - y|(l - S)) 

t2 L(l - g)2 

+ -^^y{\yHi-q)\ > S\x - y\) 
q2 



(71) 
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Let's choose 5 = exp ( — |a; — 2/|(dD(ea;-j/)\/t) ^) and q = l — 2D{ex-y)C45 
(we will use the notation D{1) := ^IDl). 

For \x — y\/\/t > ma,x{dD{ex-y)\n{4D{ex-y)C4),3dD{ex-y)) (which ba- 
sically says that the heat kernel is far from its diagonal behavior) one has 
5 < 1/10 and 1/2 < g < 1. Using the Aronson type estimate H68^ one 
controls the second term in 17H 

P.(iyt{i-,)l > 5\x-y\) < CseM-j^^f^^) (72) 

By the properties 1651 . 16611 and the corollary l2.1Ul one controls the first 
term in ||7ij: for r < 1 with r = p ^ \x - y\{l - S) - C^, and 

Ci = (2e(/i - D{cy))to)yD{e^-y) one has 

2 

P.(j/t,-e.-. >\x-y\{l- S))) < ei"-' exp ( - (1 - r^)^^^£—^) (73) 

Combining JHl and using the value of q and 5 given above 

one easily obtains the estimate 1701 of theorem 13.31 under the conditions 

□ 

Now theorem 12.81 is a straightforward application of theorem 13.31 and 
a trivial adaptation of the constants appearing in theorem 13.31 Con- 
sider p{t, X, y) the heat kernel associated to the Dirichlet form 13211 . Since 
p{t,x,y) is continuous in L°°(T'^) norm with respect to U (we refer to 
| |Chen et al., 1998 ) whose result can easily be adapted to our case) and 
C°°(T ) is dense in L°°(T'') with respect to that norm, one can assume 
U to be smooth and the general result follows by observing that the esti- 
mates in theorem 13.31 depend only on Osc(f/). 

By definition yt has symmetric probability densities with respect to the 
measure mu and the following Aronson type upper bound is available 
USeignourel, 1998| l. 

P(t,:.,y)e~^^<^' <Ce(^+'''°-(^'iexp(-^-^). (74) 

t2 4t 

It follows that the conditions 16 711 and 1681 1 are satisfied with constants 
C2, C3, C4 depending only on d and Osc(l7). Now write xi the solution of 
the associated cell problem: for I £ S , LuXi ~ —l^U with x(0) = 0. 
Using the theorem 5.4, chapter 5 of jStampacchia, 1965J on ellipt ic equa- 
tions with discontinuous coefficients (see also" ^^tampaccnia, 1966| l), using 
the periodicity of x and observing that Xi{^) = ^-^ ~ Fi{x) where Fi is 
harmonic with respect to Lu one easily obtains that 

Cx = WxAU < Cdexp((3d+2)Osc(t/)) (75) 

From Ito formula one has l.yt = x + Xiivt) ^ Xi(^) + /o(^ ^ \7xi)d'^s, 
which corresponds to the decomposition given in 16511 . The martingale 
can be written l.Mt = Jq{1 — VxO'^'^s and its bracket is equal to < 
l.M, l- M >t= J^\l -\'xiiys)\^ds. It is easy to obtain from the theorem 
3.9 of HCilbarg and Trudinger, 1983| l that 

fi = llVxdIoc < Cd{l + ||V;7|U) exp ((3d + 2) Osc(t/)) < 00 (76) 
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Writing 0; the periodic solution of the ergodicity problem Lu = I'Vxil^ — 
HD{U)l {(piiO) = 0) and observing that = - HD{U)li>i where 
Lu'^t'i = 1 it is easy to obtain from 1751 . the theorem 5.4, chapter 5 
of | |Stampacchia, 1965| l and the periodicity of 4>i that 

C0 = Moo < Cdexp((9d + 4)Osc(l7)) (77) 

Since, from the Ito formula 

E,[< l.M, l.M >t] = ^[<l>[yt) - <l>{x)\ + tHD{U)l (78) 

the martingale satisfies the conditions of theorem 13 . 31 with /2 = ^ID{U)1 
and to = C^/lyf-i — Amm(I5)). Now one can use theorem 13.31 to obtain a 
quantitative control on the heat kernel. It is important to note that all 
the constants appearing in that theorem only depend and d and Osc(C/) 

except may be fei = 30(e(/i — Amin(I)))to) ^/Amin(£^) in which /i appears. 
This is where the trick operates, indeed (/i — \min{D))to — which is a 
constant depending only on Osc(i7) and d. Thus in reality all the constants 
only depends on the dimension and on Osc(?7). Which proves the upper 
bound in theorem 12.81 

3.1.3 Lower bound estimate ()40|) of theorem 12.81 

Let yt the difi'usion associated to the Dirichlet form 13211 . As it has been 
done in subsection 13 . 1 . 21 one can prove the estimate 141)1 assuming that U 
is smooth and the general case will follow by the continuity of the heat 
kernel with respect to U in L°°{T^) norm. 
First, we will need the following estimate. 

Proposition 3.4. For I G A > fc5,ti,Osc((7) md fe6,d,Osc{!7) A < t one 
has 

nvY-l > A] > -±= / e'^ '^dz (79) 
4V27r Jx 

Proof. For I G let Fi,xi^<i>i be the functions introduced in 13. 1.21 Write 
Tt the filtration associated to Brownian motion appearing in the SDE 
solved by yt- Fi(yt) is a (P, jrt)-continuous local martingale vanishing at 
such that (Ito calculus) 

<Fi,Fi>t^tD{V) + 4>i{yt) + Mt (80) 

with Mt = — f^\7(j)i{ys)dujs- Since < Fi,Fi >oc= oo a.s. by Dambis, 
Dubins-Schwarz representation theorem Bt — Fi {yxt ) is a {J-Tt )-Brownian 
motion with Fi{yt) — B<Fi.Fi>t and 

Tt =inf{s ■.<Fi,Fi>s>t} (81) 

The idea of the proof is then to show that probability of yt to move 
away from behaves like the probability of a BM of variance D{1) to 
move away, to achieve this it will be sufficient to show that Mt becomes 
negligible in front of tD{l) using the corollarv l2.10l to control P(Mt > x). 
More precisely we will use the following lemma. 
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Lemma 3.5. For 



A>0, v>\\(t)i\\oo, /i>0, \ + \\xi\\oo + ^i < C2ii^jDiJ)tu^ (82) 
one has 

V[yt.l > A] > l/2P[Bo{i)t > A + llxdloo + m] - lP[|A/t| > - ||0i||oo] (83) 

Proof. Let A > 0, from the representation theorem P[-Fi(?/t) > A] = 
lP[^D(!)i + i^t > A] with iJt = B^Fi,Fi>t — Bu(i)f It follows that for 
M > 

P[fi(yt) > A] > P[BB(ot > A + ^i] - P[|£f I > H (84) 

It follows from ^ that for u > Q, V[\Et\ > H < ^\<t>{yt) + Mt\ > 
+ P[sup|^|^^ \BD(i)t+z - Bn(i)t\ > ^J■]■ From which one deduces 

n\Et\ > /i] < 1P[|M1 > - ||<^!||oo] + 2P[|B.| > H (85) 

Combining 1841 and 1851 1 one obtains that u > ||<?!>;||oo 

nvt.l > A] >¥[BD(i)t > A + llxi lloo + H - 4P[BD(ot 



>/i\/^^l-P[iAft| > i^-||0i||oo] 

Which leads to 18311 under the last condition in 1821 . □ 

Now let us show that 
Lemma 3.6. For Cmx < t one has 

P(Mt > a;) < 3exp(- — ) (86) 

721 

where f2 and Cm depend only on d and Osc(l7). 
Proof. Write G{x) = - ||</>i||oo0; Since 

LuG{x) = |V(^,1' - iUiWo. - 0O(|VFil' - D{1)) 
one obtains from Ito formula that 

E[<M,M>t] < 2\\(t>i\\o.E[ [ \VFi\^{ys)ds + D{l)t] + \\G\\o. 

Jo 

<2\\4,i\\^{\\(t,i\\^ + 2D{l)t)+2\\<t>i\\l,. 



(87) 



Thus Mt satisfies the conditions of the corollarv l2.10l with /2 — 4\\(l}i\\ooD{l), 
fi = |V(?!)!|^ and to = 4||(j!)i ||^/(/i — /2), which leads to 18611 by observing 
that ((/i — /2)to) 2 //2 is upper bounded by a constant depending only on 
Osc([/) and d. ' ' □ 

It follows from the equation H83|l that under the additional conditions, 
Cm{v-\\<Pi\\oo) <t, and A + ||xi||oo + A* < ^3(2^ - ||<^i||oo) (88) 
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(where C3 depends only on d and Osc((7)) one has 

Avt-l > A] > l/4P[BB(,)t > A + llx! lloo + (89) 
Choosing v = ||<^!||oo + 2/C3 (A + ||Xi||oo + m) and 



4(A+||xi||oo)5(C2\/i3(0C3t)"' 



for A > llxilloo Eind ^ > C4(d, Osc(f/))A the conditions II82II and 18811 are 
satisfied and 

M<C5(d,Osc((7))Ay|< A (90) 
and it follows from H89^ that 

nVfl > A] > \nBDii)t > A(l + Csy^) + llxilloc]. (91) 

Which proves proposition 13.41 □ 

Now, let t > 0, x,y ^ M.'' and p{t,x,y) be the heat kernel associated 
to the Dirichlet form H,S2^ . Using the chain rule one obtains that for 
< g < 1 and S > 

p{t,x,y)>Ca.Osc(u)P4yt,eB{y,5Vt)) inf p{{l - q)t, z,y) (92) 

zGB{y,SVi) 



It follows by Aronson estimates that 

p{t, X, y) >Cd,o.c{c/)lPa; {ytq G B{y, SVt)) {t{l - q)) ' 
exp ( - Cd,osc{(7),2(5^/ (1 - q)) ■ 

Now for / £ R'' let us define the probability measure Pa; as 

J-yt 



(93) 



(94) 



From now we can assume 2; := and we will fix 

l:=D{U)-'y/iqt) (95) 

and assume 

|/1 < 1. (96) 
Writing Ej; the expectation associated to ¥x one has 

Po(yt, G B{y,SVt)) =Io[e-'-^'n^^^,s(^,,^jEo[e'-^'1 

>g-yD(ur^v/{qt)-c^^osa(u),3\vls/(ithf^\^y^^ g B (y , 5^^)] Eq [e' ] . 

(97) 

Now it is trivial to check that the generator of yt with respect to ¥,x is 
L = A/2-V[7V + Z.V (98) 
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Let us write p the heat kernel associated to that generator, it is trivial 
to obtain from 19511 . Il9()l l and theorem 1.4 of ijN orris, "19971 that for z £ 
B{y,5y/i) one has 

p{tq,0,z) > Cd,osc(u)Ail^y^ exp(-Cd,osc{c/),5<5^/g) (99) 
It follows that 

^o^ytq € B{y,SVi)^ > Ca,Osc{u),aS'^q~'^^^ exp{-Cd,Osc(u),53^ /q) (100) 
Moreover for A > 

Eo[e ■^"] > mi.ytg > \]e\ (101) 

And choosing A = ID"^ {U)ltq one easily obtains from proposition l3.4l that 
there exists constants Ci, C2, C3 depending on d and Osc(?7) such that for 
\y\ > C\ and C2IJ/I < tq one has 



Eo[e'-^*1>exp(^i^^A^(l-F)) (102) 

with 

F-C3(gt/y' + |y|/(gt)) (103) 

Now let us choose 

q ■- 1 - e^Y>{-\x - y\r ) (104) 

and 

5 — (l-g)5. (105) 

With these values for q and 5 and combining Hl()2|l with H97|l and H93|l one 
obtains that for |j/ — a:| > C-j^d,Osc(u) a-nd Cg,d.osc{c7) iy ^ x\ <t one has 

Pt{x,y) > C9,d,Osc(c/)i"''''^exp(-(l - F2)\x - ylo-i(c/)/(2t)) (106) 

with 

F2 := Cio,d,OBc(t/) {t/\x - yf + \y- x\/{t)) . (107) 

It is then easy to deduce the lower bound of theorem l2.8l bv an appropriate 
shift of the constants. 



3.2 An analytical inequality for sub-harmonic func- 
tions 

3.2.1 The inequality: Theorem IZTTl 

There is no loss of generality by assuming Q, to be the segment (0, 1). We 
will give a geometrical proof theorem 12.111 explaining why we expect the 
existence of an homotopy invariant constant Cd,{i in conjecture 12. 121 The 
theorem 12.111 is proven if the inequality 1471 is true when <j) and ip are 
Green functions Gx{x,z) of — V(AV) with Dirichlet condition on d{Q, 1). 
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Let {x,y) G (0,1)^, x < y. Write fli = {z € Q : V zG{x, z)\7 zG{y, z) < 
0}. The inequahty 1471 is true if 

- [ V,G{x,z)X(z)V^G{y,z)dz < [ V ,G(x, z)X{z)V ,G(y, z) dz 

(108) 

Write {z e n : G{x,z) > G{x,y)} and Ay ^ {z G ^ : G{y,z) > 

G{x,y)}. Integrating by parts one obtains 

/ V,G{x,z)X{z)VM{y,z)dz = 0= I V,G{x,z)X(z)V,G(y,z)dz. 

J Ay 

(109) 

Now the one dimensional specificity shall be used. Since G{x, z) is increas- 
ing from to a; and decreasing from x to 1, it follows that Q.\ = {x,y) and 
[A^/ni] n {Ay/Q.i) = 0. Combining this with ([TH^ one obtains lHHSt . 
which proves the theorem. Let's note that a simple computation shows 
that the constant 3 is sharp. 



3.2.2 Equivalence with the stability of Green functions: 
Proposition 12.131 

Write fore e [0,1] A,(a:) = gU{,x)+^T[x) ^ Write V;^, the solution of -V(A,VV'e) = 

g with Dirichlet condition on and g G C°°(n), g > Q- 

Assume conjecture 12. 12l to be true, then proposition 12. 13l is proven if 

g-C,,n||T||^ 1,^1,^ <gC.,n||T||^^ (^^0) 

One obtains by differentiation (writing = — VA^V) Lx^d^tpe = —Lg^x^ipt- 
Which leads by integration by parts to 

= - / VyGxA^,y)X,{y)VyGxAy,z)T{y)g.{z)dydz. (Ill) 
Jn 

Using coniecture l2.12l 

\d.^PA < ||T||oo / \\/yGxAx,y)K{y)\/yGxAy,z)M^)dydz 

Jn (ll'^) 

< \\T\\r>cCd,ni>E- 

And integrating d^ln^pe < HrHooCti^n one obtains the upper bound in 
moi l (the lower bound being proven in a similar way). 
Conversely if conjecture 12. 121 is false one can find S > x,z G O,^ and g 
being a smooth approximation of a Dirac around z such that if 

T{y)^~Sien(yyGxAx,y)Hy)^yGxAy,^)) one has 

dAnMx) > \\T\\oo{l + S)Gd,n. (113) 
Which leads to a contradiction with 1491 . 
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3.3 Sub-difFusive behavior 

3.3.1 Exit times: theorems 12 . 11 l2T6l and proposition 12 . 71 

For r > 1, write the number of effective scales 

ne/(r) = sup{n > : i?„ < r}. (114) 

First let us prove that the exit time from 5(0, r) is controlled by the 
homogenization on those first n^f^r) scales: 

Lemma 3.7. 

—^Jr < Eo[r(0,r)] < a (115) 

D{Vo ') U D{Vq ' ) 

Proof. Write E'^, the expectation with respect to the law of probability 
associated to the generator 1/2A — Vt/V. By theorem 12.111 and proposi- 
tion 

one obtains that 

,-«o.c.{v„=:^<.Hi) < Eo[r(0,.)]/Er"'"^'[r(0,r)] < /°-(-"^.<^)+i). 

(116) 

Bounding, Un^f+i{x) by Osc(l7„) < Kq and for k > n^f + 2, Uk{x) by 
II V?7fe||oo|a;| < Ki\x\/Rk one obtains that for x G 5(0, r) 

|V;^,{r) + l(a;)| < ^0 + i^l/(p.nin " 1). (117) 

Writing p^f corresponds to the maximum number of periods of the scale 
n^f included in the segment [0,r]: Pes{r) — sup{p > 1 : pRn^f(r) < 
one obtains 

Eo° [r(0,p./(r)7?„^^.„)] 

< [r(0,r)] < E^," [r(0, (pe/(r) + l)i?„,,(.))]. 

Using E^°°'"[T(o,fcfi ( ))] = (fci?„^^M)7i3(\/o"=^<''>), dna, era 

and 11181 1 one obtains 11151 1. □ 

We will need the following mixing lemma 
Lemma 3.8. Let (gj) G {C\T^f and R e N* 

g{x)f(Rx)dx- [ g{x)dx [ f{x)dx\ < \\Vg\\^/R [ \f\dx 

Proof. The proof follows trivially from the following equation 

g{x)f{Rx)dx— / g{x)dx / f{x)dx 



f{Rx + y){g{x + y/R)-g{x)). 

j/e[0,l]<*,a;eT'' 

(119) 

□ 
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From lemma 13.811 we will deduce a quantitative estimate on the multi- 
scale effective diffusivities: 

Lemma 3.9. 

(A„,i„e-^^'i/^""")" < D{V--') < (A„,axe*^i/^-")". (f20) 

Proof. The proof of 112UII is based on the following functional mixing 
estimate (obtained from lemma 1^^ : for U,W € C^(T) and R € N* one 
has 

g-||VW||oo/-R <- 

(121) 



Then by the explicit formula I20II and a straightforward induction on n 
one obtains that (using Q) 



n-l 



fe = 



n (e'*^i'''-'= / e'^^(">da; / e^'^^^'^da;) < D(Vo"~') (122) 



n — 1 „ 

D{Vr^) < n (e-*^i/'-'= / e^^^f^'di / e-2^^(^)dx)-\ (123) 

Which leads to (IHITI by © and □ 

Combining 11201 with II115II . 111411 and lHUIl . one obtains theorem ITTl 
When the medium is self-similar, we will need the following lemma 

Lemma 3.10. 

lim -- In (DiV-^)) = Vp[2U) + Vp{-2U). (124) 

n— »oo n 

Proof. The limit (11241 is a direct consequence of the following theorem 
that is an application of the theory of level-3 large deviations (we refer to 
( |Ellis, 1985} for a sufficient reminder). 

Theorem 3.11. Let U £ C"(T'') (Holder continuous with exponent a > 
0). Let ReN, R> 2. Then 

/■ ""^ 

lim -In / expiJ" U{R''x))dx ^Pr{U). (125) 
n Jjd ' ^ 

We have written Vr is the pressure associated to the scaling shift 
induced by R on the torus: For R G N/{0, 1} one can see the torus as a 
shift space equipped with the transformation sr 

SR-.T'^ — > T'' 

- „ ^ a^'^+i (126) 

fc=i fe=i 



where for each k, x^ is a vector in B = {0, 1, . . . , ii — 1}'* and for each i G 

R*" 



{1, . . . , d} X^fcLi "Rk is the expression of Xi in base i? (i*^ G {0, . . . , i?— 1}). 
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Give B with the discrete topology and B with the product topology. 
Write fi the probability measure on B affecting identical weight l/R'^ to 
each element of B and write the associated product measure on B*' . 
With respect to the probability space (B'* ,B{B^'^ ^he coordinate 

representation process x = (a;^, . . . , , . . .) is a sequence of i.i.d. random 
variables distributed by p. When x is seen as an element of the torus T"^ 
then the probability measure induced by ^ on the torus is the Lebesgue 
measure. 

Define the empirical measure _E„ associated to the process x by 

^ n — 1 

^"(^'•) = -E'^-i^cycM.,™) (127) 

fc=0 

where cycle(a::,n) is the periodic point in B'*' obtained by repeating 
(x^ , . . . ,a;") periodically. For each i, En{x, .) is an element of the space 
A4{B^ ) of measures on B^ and invariant by the shift sr. 
Then by theorem 9.1.1 of ( |Ellis, 1985^ , {Q'^^}, the distribution on 
M{B^ ) of the empirical process {En} has a large deviation property 
with speed n and entropy function 7^^ . 

We remind that for P G M^B^'), = Jgm I^J''' {P)dP where P 

1 (2) 

denotes the marginal distribution of x associated to P and is the 
relative entropy of P with respect to fi: (rf) = In d^. 
Choosing U G C(T'^), Holder continuous with exponent a, one deduces 
from the large deviation property of {QI^'} and Varadhan's theorem that 

limiln/ eyipinE.a{x,U))dx = VR(U) (128) 

n— >oo n Jjd 

where Vr{U) is the pressure of U . We remind that 

Vr{U)= sup { [ U dP - li:'\P)} (129) 

where A^s^(-B'* ) is the space of measures on B" invariant by the shift 

Sr. 

Since U is Holder continuous 

n — l n — 1 „ 

\nE„ix,U)^Y.UiR''^)\<j:(^r 

fc=o fc=o ^ ^^3Q^ 

fc=0 

And one obtains theorem l3. Ill from II13UII and 112811 □ 



Combining H124^ with l|115|l and Ijll4|l . one obtains theorem ITi 
Now let us prove proposition 12.71 The basic properties of the pressure 
can be found in ( [Keller, 1998^ theorem 4.1.10. (note that the definition 
of the pressure given here differs from the standard one of the topological 
pressure by a constant that is dln_R, here 'Pij(O) = 0). Let's remind that 
Vr is a convex function on the space of upper semi continuous functions 
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on the torus to [—00, 00) thus VRiU) + Pr{—U) > 0. 

We wiU now remind the strict convexity of the topological pressure on a 

well defined equivalence space: 

To sr is associated a scaling operator Sr acting on the periodic continuous 
functions on T'' 

Sr : c{r) ar) 

{x~.f{x))~-.{x-.f{sRx)=f{Rx)). 

Write Ts^(T'*) the closed subspace of C(T'*) generated by the elements 
V - S^V with V &C{T'') an d fc £ N. Writ e [U] the equivalence class of 
U, then by proposition 4.7 of ( |Ruelle, 1978^ the function 

rR:C{r)/isA'^'')-^[-<^,+<^) ^32^ 

[U] Pr{U) 

is well defined on the set of equivalence classes induced by Ts^(T'*) on 
C(T''). Moreover it is strictly convex on the subset 

{lU]eC{T'')/Isn{^'') ■■ f U{x)dx^O}. (133) 

We will now prove proposition 12.71 since for c G R, 1^(11 + c) = ViU) + c, 
it is sufficient to assume J^^ U{x)dx = Q and show that 

n — 1 

Vr{2U)+Vr{-2U) = ^ lim -|| V Sj^kUW^ = 0. (134) 

n — >oo 72 — ^ 

(<=): This implication is easy since 

Q<Vr{2U)+Vr{-2U)< fim -|| V Sflfct/||oo. (135) 

n — >oo n — ^ 

fe = 

(=>): Assume Vr{2U) + Vr{-2U) = then let e > 0. Then by the strict 
convexity of the pressure as described above there exists W\ , . . . , Wk G 
C(T'')andmi,...,mfc G N/{0, 1}, Ai, . . . , Afc G i? such that W = Ep^i ^p(W^p-^ 
SRm.pWp) and \\U — W^||oc < e. Since X]p=o SrpW remains bounded it 
follows that 

n — 1 

lim -II y Sfl.fJIloo < e. (136) 
which leads to the proof. 



n— >oo Tl 



3.3.2 Mean squared displacement: proposition 12.51 theo- 
rem I 



Let j/t be the solution of (Q. Write 

n;iu{t) ^sup{n£n : Rl<t} (137) 
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Ufiu shall be the number of fluctuating scales that have an influence on 
the mean squared displacement at the time t (the effective scales plus the 
perturbation scales). Chose the number of perturbation scales to be 

riper = mi{n e N : i?^^.,^^_„e""-^nO* < (138) 

We will now prove the following proposition 

Proposition 3.12. For pmin > lOe'^^'^^^ and t > Rg, Uper is well defined 
and 

Cie-^"^^^'^"D{Vo''")t < E[y?] < Cze'^"^"-'^" D{Vo'''')t. (139) 

Proof. The proof of 11391 1 is based on analytical inequalities that allow to 
control the stability of the homogenization process on the smaller scales 
under the perturbation of larger ones. More precisely we will first work 
on an abstract decomposition of V given by (|5J into effective scales U 
perturbation scales P and drift scales T: V ^U + P + T with {U, P, T) G 
C°°(T^^)xC°°(Ti^^)xC°°(R), Ru,Rw& N/{0, 1}, Rw/Ru = Rp e 
and W = U + P shall correspond to fluctuating scales. 
Write the solution of the cell problem associated to Lw {LwXw = 
— VTy, Xw(0) = 0) and Fw{x) = x — x^ {x). Since Fw is harmonic with 
respect to Lw = Lv + VTV one obtains by Ito formula that Fw(yt) = 
V Fw{ys)dLOs — /(* VTV Fw{ys)ds from which one obtains that 

^ is\ 

(140) 



{i^-t\\VT\\L)ni \VFw[ysTds\ 



< E[Fl-{yt)] < 2(1 + i||Vr||L)E[ / \VFw{ys)\^ds]. 

Jo 

Write x^ the solution of the cell problem associated to Lp and F^ = 
X — x^ ■ We will show that 

Lemma 3.13. Fw = Fp - with 



and 

< 2(1 + A\\VP\\a.)e^°''^''^Rw/Rp. (142) 

Proof. The inequality 1141^ is a direct consequence of the explicit formula 
F^{x) = Rw e^^^y^dy/ J^"^ e^^^^^dy. The inequality follows 
from the explicit formula 



H"{x) = Rv 



^ /o^* e'^Piy)dy /q^"' e'^iP'^y^+umdy 
noticing that the period of P and U are Rw and Rw/Rp and lemma 



□ 



The long time behavior of E[J^ \VFw{ys)\'^ds] is a perturbation of 
D{W)t as shown in the following lemma 
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Lemma 3.14. // 

RP > 16e*°=^(^HllVP|U + |lVT|U)e^ll^^l'~/^- 
then for t > 

+ (l/6)e~'*°="'-^^L>(W)t < E[ /' lVF'''T(y,)ds] 

Jo 

and 

E[ /* lVF'^l'(y,)ds] < 6e''°='=(-^'D(VK)t 
Jo 



(143) 



(144) 



Proof. For the proof of I I143II and 114411 by scaling one can assume that 
Rw ^1 and Ru = 1/Rp. Write for C > 

rx 2V(y) ry 2(P-T){z) ry -2{P+T){z) 

(pC=2 —r^ / dz-C dz\dy. 

Jo e^^'iy^dyUo J^e^P(-)dz Jo e-^^^dz J 

(145) 

Using lemma to separates the scales in 111451 . it is an easy exercise to 
obtain that if Rp > IBe-^^^'f^'dl VP||oo + |1 Vr||oo)e''ll^^ll°°/-^^ then 

. for C = 6e40=^(^) one has sup. < 900 ^"°;°"'' e""^^"-/^ 
. for C = -"'7'"' one has inf m 0c > -lOO^i^^^e^ll^^"-/^ 

Observing that iv<^c = K - - CDiW) one deduces llT^ and llT^ 

by applying Ito formula. □ 

Combining ifTIH)! . ifTlTji . itTl^ . ifHi^ and choosing [/ = ^^"/i" ^ 

and nyju as defined in 113711 one obtains that for pmin > Cki.Kq 

D(l/o"^'")te-«"-'-^V24 - i?^^,^_„^^^500e«"--^« < E[y?], (146) 



nvt] < {DK^'^t + ^n.^^-n.^Je^-^-^^SOO. (147) 
Which leads to 11391 by the choice 11381 for riper- □ 

By the uniform control of the ratios pOjl one obtains quantitative 
estimates on the number of fluctuating and perturbation scales 11371 and 
11381 : combining them with the control 11391 and the exponential speed 
of convergence of the multi-scale effective diffusivities towards zero lfT20|l . 
one obtains proposition 12.51 and theorem 12.21 
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3.3.3 Heat kernel tail: theorem !^ 



As it has been done for the mean squared displacement, the proof of 
theorem 12 . 31 shall follow from an abstract decomposition of the potential 
V. More precisely, let Rw G N/{0, 1}, (W,r) G C°°(TA^^) x C°°(K), 
(II VT||oo < oo) and write V = W+T and yt the diffusion associated to Ly- 
It has been shown in the proof of proposition 13 . 1 2l that by decomposing 
W into U + P where U is of period Rw / Rp G N one has for all t > and 
aU a; G R'' 

E.[ |VF^|2(y.)ds] < C2D(W)t + (148) 
Jo ^p 

where the constants Cj , (^2 are those given by the equation 11441 . We will 
now show that from the control 114811 (and Hx'^Hcx) < Rw that is given 
by the explicit formula of the solution of the cell problem) one can deduce 
the following lemma: 

Lemma 3.15. For 

Rw < h/2 (149) 



||VT||oo2'(C2i:>(W^)) ^ < {h/t) < {Rp/{Rw\[ci))(,2D{W) (150) 

and 

I— h2 

{Rp / {Rw\l Ct))C,2D{W)e ^''<20(w)t < (/j/^) (151) 

one has 

P[yt >h]<Ce 2'JC2-D(w)t . (152) 

Proof. The proof of 11521 1 is based on a control of the Laplace transform 
of yt, more precisely it is well known that for A > 0, and h > one has 
V[yt >h]< E[e^(^*-''>]. Observing that yt = iVt) + Jo VF^{jj,)dLUs - 
J^VT.VF^{ys)ds and using \\x ||oo ^ Flw one deduces by the Cauchy 
Schwartz inequality that 

,1 (153) 

]gjg2v^||VT||„oA(/o* |VF,»'(a,)|2ds) 

If X is a positive bounded random variable, ^' > and A' > it is easy 
to show by integrating by part over dF{X > x) and using W{X > x) < 
E[exp(A'(X - x))] that 

E[exp(M'^/X)] < l + ^'exp(^^)yjE[exp(A'X)] 

Applying this inequality to IT^ with X = [VFi^"^' (ys)fds, A' = 8A^ 
and/i' = 2Ayi||VT||^ and observing by Ito formula that E[e2^JoV^r(!/e)<^'^=] < 

>h]< Ce^{«w-Mgl|VT||Lt/4^j^8AVo'|VF»'(,.)|^dsj 
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Now observe that V (ys)dLUs satisfies the conditions of theorem 12.91 



with /a = C2D(W), and to(/i - /2) = -I^C?. It foUows that for 

SX" < {R%) / {2eR'wCt) (154) 

one has 

Assuming Rw < h/2 and choosing A = 32(^3D(vy)t the condition on A in 
11541 1 is satisfied under the right inequality in 115011 and one obtains 

P[yt >h]< Ce'^'i20(w)t^\\^Tf^i/'i{^^ip(^c,2D{W)if) /{h^iCtfRw)- 

From this the result l|152|l follows easily by assuming the left inequality in 
11501 1 (that basically says that the influence of the drift scales ||Vr||oo is 
small in front of the influence of the fluctuating scales) and the condition 

iirni . □ 

Now let's choose W = V^''-, P = V;7;;„„^^„+i, T = (Rw = 

Rnfi^, Rp ~ Rnji^l Rnfi^-nj,„) in lemma 1^.151 For p £ W define the 
function 

np..(p) = inf{n G N : [R^ / Rp.„)e-^'"'° D{Vr^)'^ > 2" (155) 

nper (p) corresponds to the number of perturbation scales among p fluctu- 
ating scales. We will from now assume that pmin > 2^e^^^i , which implies 
that Uper is well deflned and 1 < Uper (p) < P- Define 

Ufiuit/h) = inf{n G N : 2«(Xi/i?„+i)e2"'"='-(")-^« (L»(l/o"))* < h/t} 

(156) 

nfiu — nper corresponds to the number of fully homogenized scales given 
t/h. Ufiu is well defined and greater than 1 under the following assumption 
that basically says that homogenization has started on at least the first 
scale. 

{R2/K^)e^'^''2-'' <t/h (157) 

By the definition of nfi^ the left inequality in 115011 is satisfied. Using 
115511 . the right inequality in 11501 is implied by the definition of n^er- 
The inequality 11491 is satisfied if 2i?„^.,^^ < h; by the definition of Ufiu 
this is implied by the following inequality that basically says that the heat 
kernel behavior is far from its diagonal regime. 

/{D{Vo^'^)^t) > 2A'ie''-^"2''e''"''='-("^'"^^° (158) 

By the definition of nfiu and Uper, the inequality 115111 is satisfied by 
the following inequality that also says that the heat kernel is far from its 
diagonal regime. 



2l4g4(nper + l)ifo 



[i?n,,„+l] <ft7(m"''"W (159) 
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With this assumption, it follows by the inequality 11521 that 



¥[yt>h]<Ce 2ii.^"p-'^o,,(v;/'")t (^gg) 

Using the control iHUli on DiVo^'"), and (O on the ratios one obtains 
theorem l2.3l The condition 11571 is translated into the first inequality in 
1161 and the conditions 11581 . 11591 into the second one. 
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